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Abstract 



We study resonant Andreev tunneling through a strongly interacting quan- 
tum dot connected to a normal and to a superconducting lead. We obtain 
a formula for the Andreev current and apply it to discuss the linear and 
non-linear transport in the nonperturbative regime, where the effects of the 
Kondo resonance on the two particle tunneling arise. In particular we notice 
an enhancement of the Kondo anomaly in the I — V characteristics due to the 
superconducting electrode. 
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The way a superconducting electrode affects the electrical resistance of a small meso- 
scopic region via Andreev tunneling, has been the subject of an impressive research activity 
over the last years This has lead to the discovery of a number of new phenomena. 
Examples are the zero-bias anomalies in Normal-Insulating-Superconducting (NIS) struc- 
tures or temperature (and voltage) re-entrant behaviour in the conductance of a small 
metallic wire in contact with a superconducting lead 0. Several of these phenomena have 
been successfully explained by models of non-interacting electrons Electron-electron 
interaction strongly modifies the two-particle tunneling. In small capacitance junctions, for 
example, it leads to the Coulomb blockade of Andreev tunneling |TI]], whereas in a Lut- 



tinger liquid-Superconducting structure anomalous I — V characteristics appear due to the 



excitation of the low lying modes ||Tl 



Electrical transport via resonant tunneling in strongly interacting systems has been in- 
tensively investigated because of the possibility to study the behavior of an Anderson-like 
impurity model. The prototype model is a Quantum Dot (QD), which mimics the impurity 
level, coupled by tunnel junctions to two non-interacting leads. When the dot level is far be- 
low the Fermi energy of the leads, the formation of a spin singlet between the impurity spin 
and the conduction electrons gives rise to a many-body resonance at the Fermi energy. Such 
a resonance, known as Abrikosov-Suhl or Kondo resonance, manifests in a peak in the den- 



sity of states at the Fermi level and leads to a perfect transmission at zero temperature |]T2 
The two reservoirs, differently from the case of the Anderson impurity in metals, can be kept 
at different chemical potentials. Nonlinear transport through a QD provides then a beatiful 



tool to study a Anderson impurity out of equilibrium [|T3[. The nonlinear I — V characteris 



tics and the Kondo resonance in nonequilibrium conditions have been studied by a variety 



of methods: the equation of motion for the Green's function |T3], variational procedure [14 



a newly developed diagrammatic technique and the 1/N expansion |TB|. Experimental 



evidences of Kondo-like correlations have been reported in metal point contact [T^ and very 
recently in SET transistors |]TB|. Zero bias anomalies have also been reported recently in 
experiments involving superconductor/Anderson-insulator interfaces, though their physical 
interpretation is still under debate ||T9| . 

Remarkably a QD offers the, yet unexplored, possibility to study a quantum impurity 
coupled to two different types of electronic systems, such as superconductors and normal 
metals at the same time. The development of superconductor-semiconductor integration 
technology should allow to study the signatures of the Kondo effect in the Andreev tunneling. 

In this Letter we will study Resonant Andreev Tunneling (RAT) in a QD coupled to a 
normal and to a superconducting reservoirs as shown schematically in Fig. |l|. Differently 
from the conventional problem of a Kondo impurity in a superconductor here the low 
lying excitations are still present due to the presence of the normal electrode. 

The model Hamiltonian for the N-QD-S system under consideration can be written as 

H = Hpf + Hs + Hd + Ht^n + Ht^s (1) 

where Hn = Ek,a ^kCN,k,aCN,k,'^ , Hs = J2k,a ^kCs,k,aCs,k,^ + Ek{^cl^k,^cl^-k,i + ^■^■) and 
Hd = ^ddlda + f/?7.d,|nrf,j are the Hamiltonians of the normal lead, the superconducting lead 
(A is the superconducting gap) and the dot respectively. The single particle energy is 
double degenerate in the spin index a and the interaction is included through the on-site 
repulsion U. The position of the dot level can be modulated by an external gate voltage. 
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The tunnehng between the leads and the dot is described by Hx^n = Y.k,cr{Vv,k,(T'^lj,k,a,d'CT + c.c.) 
where r] = N, S and K?,fe,cr is the tunnehng amphtude. 

The average current, which for convenience we compute in the normal electrode, is given 
by J = 2eIniJ2k < '^N,k^zHk^ >, where we have adopted the Nambu notation (the hat 
indicates matrices in the Nambu space): '^N,k = {cN,k,h^N,-k,i) $ = {d^,d\). Hk is a 
diagonal matrix with elements Hn = VAr,A;,T,-f^22 = V^_j.^^. 



By means of the Keldysh technique, as employed in Ref. pT[, the current / can be 
rewritten in the form 

-r^Tr{f,G^(e)[S«(e)/^(e) - /^(e)S^(e) + S<(e)]G^(e)} (2) 

-OO ZTT 

where G^^^^ , G"^ are the retarded (advanced) and the lesser Green's functions of the dot (for 
example G^{t) = -iOit) < {$(t), $t(0)} >). In deriving eq.(D, the relation G< = G^t<G^ 
has been used. Here r7v(e) = 2TrJ2k |^A',fc,crp5(e — ^k)- For the sake of simphcity, from now 
on we will neglect the spin and energy dependence of the tunneling matrix elements. The 
diagonal matrix has elements fN,ii = /((e + eV^)/T) and fN,22 = 1 — + eV)/T) 

if the normal electrode is kept at a chemical potential which differs by eV from that of the 
superconductor (/(x) is the Fermi function and V the voltage drop) . 

To obtain a suitable expression for the current, we need to evaluate the lesser self-energy. 
In the non-interacting case, S^(e) can be computed exactly and it is expressed in terms of 
the retarded and advanced self-energies as Sq (e) = —J2ri=N,s(X'o,r]i^)fr]{^) ~ fvi^)^o,v(^)) ■ 
In the interacting case, we generalize Ng's ansatz to the present case [^. The lesser and 
greater self-energies are assumed to be of the form 

S< = , S> = ±>A (3) 

where A is a matrix to be determined by the condition S< — S> = S-^ — S"^. This ansatz 
is exact both in the non-interacting limit, U = 0, and in the absence of superconductivity, 
A = 0. Moreover it guarantees automatically a current conserving scheme for any given 
approximation procedure to evaluate the retarded Green's function. As a result we obtain 
= Sq (S^ — S^)~^(S^ — T,^) . The expression for the current can be greatly simplified 
in the relevant limit U,A^ ksT, V. In this case the non- interacting self-energy due to the 
superconducting lead Sq s^'' is real and purely off-diagonal whereas that due to the normal 

lead, S^Jv^ is diagonal. As a consequence, we obtain the following form for the Andreev 
current through a QD 

-r^Tr{f,G^(e) S«(e), Me)] G%)} . (4) 

-OO ZTT L J 

This is the first central result of this work. It generalizes to the case of a strongly interacting 
dot the formula valid in the non interacting case and allows to study the transport 
through a N-QD-S under non-equilibrium situations and in the nonperturbative regime. 

The presence of a superconducting lead introduces, besides U, another energy scale in 
the problem, the energy gap A. One is often interested to voltages and temperatures well 
below the energy scale set by A and could be tempted to take the A ^ oo from the outset. 
In an interacting quantum dot, the A ^ oo limit cannot be taken, in contrast to the non- 
interacting case. In order to have coherent Andreev scattering, the two electrons enter the 
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superconductor without double occupying the QD (which is forbidden in the U ^ oc hmit). 
This can happen only on a time scale of the order 1/A. We will then consider the case 

The generalization of the decoupling scheme for in the presence of superconductivity 
is conveniently done by rewriting the superconducting lead Hamiltonian, Hs, in terms of 
quasiparticles operators by means of a Bogoliubov transformation. The reason for trans- 
forming to the quasiparticles basis is dictated by the type of approximations introduced in 
the equation of motion approach. In fact, one may generalize to the case of a supercon- 
ducting lead, the decouphng usually introduced in the normal case. Such decoupling, which 
neglects correlations in the lead, is done in terms of the equilibrium number of quasiparticles. 
The final expression for the matrix Green's function of the dot is 



5,26 
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(5) 



where (Tat h is the U = oo-limit self-energy considered in Ref. 



e — Cfc + iO' 



(6) 



and (7s 11 is the corresponding self-energy due to the superconducting electrode 



EM 



e + Ek + iO+ 



(7) 



{vl = 1/2(1 — ek/Ek) and Ek = \j t\ + |AP). The self-energy for the hole propagator can be 
obtained using the property cr7v(5),22(e) = ^*N{s),iii~^)- '^^^ self-energy as,ii is weakly energy 
dependent due to the low energy cutoff provided by A (in what follows, we consider uic 3> 
A ^ T, e with uo^ the bandwidth). At energies much smaller than the gap quasiparticles 
present in the dot cannot decay by tunneling into the superconductor (as the imaginary 
part of the diagonal elements of the self-energy as vanishes). As a result the contribution of 
the self energy due to the superconducting lead as,ii simply shifts the dot level to the new 
value Q ~ erf + (r5/27r) IncUc/A. The divergence of the level energy renormalization with 
A reveals that the process occurs via a virtual state in which a quasiparticle is created in 
the superconductor. It is however due to the weak (logarithmic) nature of this divergence 
that RAT can can be observed in any realistic situation. Note also that the off-diagonal 
element, T s/2, entering eq.(|^), has the same form of the non interacting case apart from an 
overall phase-shift of vr. This vr phase-shift comes from a relative cancellation between the 
non-interacting and interacting off-diagonal self-energies. Substituting the espression of the 
QD's Green's function in Eq.(|) we get the desidered result 



de 



f(e-eV)-f{e + eV) 
2e 



G 



(8) 



with 
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h [4(e - ei)(e + e^) - ri,^r2,,v - r|]' + 4 [ri,;v(e + £2) + r2,^(e - ei)]' 

where e^2) = + ReaN,ni22) + Reas,ii(22), ^i{2),n = -2 ImaN, 11(22) ■ Eqs.(|J|) are the main 
results of the present paper. 

The spectral function Gns{^) associated with the resonant Andreev tunneling is plotted 
in Fig.0 for various bias voltages (for comparison, the non-interacting case is shown in the 
inset). Several features are worth noticing. First, two peaks at ±e are due to particle 
and hole bare levels. Note that in the interacting case the bare level energy includes the 
renormalization due to the superconducting electrode self-energy as discussed above. Second, 
at low temperatures a Kondo peak develops at the Fermi energy. Quite remarkably, at finite 
positive (negative) voltages the Kondo peak shifts pinned to the Fermi level of the normal 
metal while a small kink develops at negative (positive) voltages. At finite voltages hole 
and particle energies differ by 2eV, and while the electron (hole) is on resonance for positive 
(negative) voltage, the Andreev reflected hole (electron) is off resonance with respect to 
the shifted Fermi level. Compared to the N-QD-N case, here there is no peak splitting. 
Third, the Kondo peak remains rather pronounced even in the nonequilibrium situation. 
The differences between the N-QD-N and the N-QD-S cases, can be traced back to the 
fact that the superconducting electrode acts simply as a boundary condition, even in the 
nonequilibrium situation, the Kondo resonance being achieved through the tunneling into 
the normal electrode. For this reason the Kondo peak is shifted but not suppressed. 

The differential conductance, for various temperatures, is shown in Fig. ^. In the low 
temperature regime the zero bias anomaly associated with the Kondo effect is clearly seen. 
We note that the Kondo peak survives up to temperatures of the order of r/4, about one 
order of magnitude larger than the N-QD-N case. This effects is shown in the inset of 
Fig. ^ where the zero bias anomaly for the N-QD-S (present work) is compared with the 
corresponding quantity in the N-QD-N, Ref. The anomaly in the RAT clearly survives 
in a temperature regime where is already absent for normal tunneling. This may be useful 
for the experimental detection of the effect. 

In this Letter we studied the RAT in a Normal metal - Quantum Dot - Superconductor 
device. An explicit form of the current through the device is obtained. The analysis of 
the I — V characteristics has been carried out in the limit U ^ A ^ T,V, where two 
electrons tunnel almost indipendently since the onsite Coulomb repulsion (which is the 
largest energy scale in the problem) prohibits the double occupancy of the Dot. In this case 
the Kondo effect enhances the Andreev conductance at low temperatures. In the opposite 
limit, A ^ U ^ T,V, different processes dominate the transport and a suppression of the 
Andreev tunneling is expected. The results were obtained by generalizing, in the presence of 
superconductivity, an ansatz due to Ng and by means of the equation of motion method to 
determine the Green's function of the Dot. While the expression for the current in Eq. (^ 
holds under general circumstances, the results for the spectral function and the differential 
conductance are quantitatively valid for temperatures larger than the Kondo temperatures 



and only qualitatively valid for lower temperatures |ll3| , ^ . We then expect that, even though 
a more refined treatment is necessary for more quantitative results in the Kondo region, all 
qualitative features of our findings will survive. The combined effect of a finite A and U is 
currently being investigated and it will be a subject of a forthcoming publication. 
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FIGURES 



FIG. 1. The system under consideration. A Quantum Dot coupled by tunnel barriers to a 
Normal and to a Superconducting electrodes. The position of the level in the dots can be tuned 
by means of the gate voltage Vg 

FIG. 2. The spectral density for the two particle tunneling G^si^) is plotted for various bias 
voltages {V = solid hue, V = ±0.005 dotted line, V = ±0.01 dashed line, V = ±0.015 dot-dashed 
line, ed = —0.07, F^ = F^v = 0.02 and T = 0.0001, in units of the bandwidth coc)- In the inset the 
non interacting case is shown for comparison. 

FIG. 3. The differential conductance of the NQDS device, in units of 4e^//i, is plotted for 
different temperatures (T = 0.0001 solid line, T = 0.001 dotted line, T = 0.01 dot-dashed line, 
€fi = —0.04, Fg = Fat = 0.02, in units of the bandwidth lOc)- In the inset the zero bias anomaly, 
detected by measuring the relative change of the differential conductance G{V) between zero and 
high voltage {AG/G = (G(0) - G(-0.04))/G(0)), is shown for the N-QD-S (triangles ) and for the 
N-QD-N (filled circles) [|l| 
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